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Abstract 

We derive the Hamiltonian associated to a quantum stochastic flow by 
extending the Albeverio-Kurasov construction of self-adjoint extensions 
to finite rank perturbations of nonsemibounded operators to Fock space. 

1 Introduction 

To analyze a dynamical system with Hamiltonian K = Kq + T, occurring as a 
perturbation of a free Hamiltonian Kq, we transfer to the interaction picture by 
means of the wave operators 

V {t) = Uo (-i) U (t) 

where U and Uq are the strongly continuous groups generated by K and Kq 
respectively. The family V = {V (t) : t > 0} then satisfies the differential 
equation 

V{t)^-^Tit)V{t) 

where T (i) = Uo{—t)TUo{t). However, V is not a semi-group, but instead 
forms a C/o-cocycle, that is ^ (< -I- s) = Uq (-t) V (s) Uq (t) V {t) for aU t,s>0. 

Conversely, given a strongly continuous group Uq and a strongly continuous 
t/g-cocycle V, then we may define a family U by 

rru^f UQ{t)V{t), t>0; 

^ ^ I V{-t)^UQ{t), t<0. 

and this is readily seen to be a strongly continuous group. Consequently we 
deduce the existence of an associated Hamiltonian K. 

The somewhat surprising feature of the converse is that V does not need 
to be strongly differentiable. In particular, it applies to the class of quantum 
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stochastic evolutions modelling open system dynamics. In such cases, the exis- 
tence of Hamiltonians Kq and K is immediately apparent, however, it must of 
necessity be the case that domains do not have dense intersection and so the 
formal subtraction T = K — Kq does not lead to a self-adjoint operator in any 
meaningful way. 

A well studied model is the dilations of an irreversible dynamics using Fock 
space over L^-functions of time where the free dynamics is second quantization 
of the time shift. It has been a long standing problem to characterize the associ- 
ated Hamiltonian K for these models [T]. The major breakthrough came in 1997 
when A.N. Chebotarev solved this problem for the class of quantum stochastic 
evolutions satisfying Hudson-Parthasarathy differential equations with bounded 
commuting system coefBcients [2] , [8] . His insight was based on scattering theory 
of a one-dimensional system with a Dirac potential, say, with formal Hamilto- 
nian 

k^id + ES 

describing a one-dimensional particle propagating along the negative x-axis with 
a delta potential of strength E at the origin. (In Chebotarev's analysis the S- 
function is approximated by a sequence of regular functions, and a strong resol- 
vent limit is performed.) The mathematical techniques used in this approach 
were subsequently generalized by Gregoratti [14] to relax the commutativity 
condition. More recently, the analysis has been further extended to treat un- 
bounded coefficients [T7j. 

Independently, several authors have been engaged in the program of describ- 
ing the Hamiltonian nature of quantum stochastic evolutions by interpreting 
the time-dependent function T {t) as being an expression involving quantum 
white noises satisfying a singular CCR [9] [10] [11] [3] . This would naturally sug- 
gest that T should be interpreted as a sesquilinear expression in these noises 
at time t = 0. Specifically one considered time dependent Wick ordered ex- 
pressions T {t) = J2ij (0 flj (t) + ^Mal (*) + Eojaj {t) + Eqo where 
the Eij = ijjj are operators on the initial space and {t) , aJ (t) are delta- 
correlated quantum white noises corresponding to the formal derivatives of the 
annihilation and creation, process. With (t) = Uq {—t) ai (0) Uq (t), we see 
that T {t) = Uo{-t)TUo{t) where T = T(0). We shall show that this in- 
tuition essentially provides the correct answer, though without using quantum 
white noise! This approach inspired W. von Waldenfels to give an alternative 
construction of the associated Hamiltonian, for diffusions 18] and simple jump 
processes |19| . however this was formulated through the conventions of kernel 
calculus. 

The aim of the current paper is to complete this program by returning 
to the original one-dimensional model considered by Chebotarev. We alterna- 
tively consider this as a problem of finding a suitable self-adjoint for the singular 
Hamiltonian [16| . Here the generator of the free dynamics fco — id is not semi- 
bounded and the J-perturbation is viewed as a singular rank-one perturbation. 
We employ methods introduced by Albeverio and Kurasov [J [S] [5] to construct 
self-adjoint extensions of such models. In particular, the boundary conditions 
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that ones imposes at the origin corresponds to a phase change of s = f , 

which should be contrasted with the condition s — e^*^ deduced by Cheb- 
otarev. We show that the constructions of Albeverio and Kurasov are amenable 
to second quantization and the form of T suggested by quantum white noise 
analysis is precisely what is needed to obtain the description of the associated 
Hamiltonian K derived by Chebotarev, Gregoratti and von Waldenfels. The 
construction avoids both the use of quantum white noises and the unwieldy 
complexity of the kernel calculus by instead using the defect vectors which lie 
in the underlying one-particle Hilbert space. One important subtly is that the 
continuous tensor product decomposition for Fock spaces is here implemented 
by the Sobolev space W^'^ inner product and not the usual inner product. 

We show that there is a natural class of boundary conditions giving rise 
to different extensions, and therefore different physical representations of the 
same problem. These are parameterized by what effectively are the additional 
complex damping terms (describing for instance energy shifts, e.g., the Lamb 
shift [2]) and which have earlier been referred to as a degree of "gauge freedom" . 



2 Singular Perturbations 

We shall adopt the standard convention that the sesquilinear Hilbert space inner 
product is conjugate linear in the first argument and linear in the second. Let ^ 
be a conjugate linear functional on some domain of functions on the real line. Its 
adjoint is the linear functional on the same domain defined through complex 
conjugation, that is, — ^ (0)*. With a standard abuse of the Dirac bra-ket 

notation, we shall write = {^\4')i(. (</>) = (010: O'' more simply 

even though the functionals need not correspond to vectors in the Hilbert space 
L2 (M). 

The weak derivative of a measurable function (j) will be denoted as d(f). For a 
fixed Sjhe a Hilbert space and / an open subset of K*^, we shall write W^''^ {I, S^) 
for the Sobolev space of i^-valucd functions (j) possessing a (weak) derivative and 
such that both such that (p and d4> are square-integrable. This is again a Hilbert 
space with inner product 

m)i,2^ |^(0*v+90*9v^). 

Note that the corresponding norm ||-||-^ 2 is then the graph norm associated with 
the first derivative operator d. 

Formally one may consider the Hamiltonian 

describing a one-dimensional particle propagating along the negative a;-axis with 
a delta potential of strength E at the origin. As a mathematical problem, we 
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then have a perturbation of a sclf-adjoint operator ko — id, which is not semi- 
bounded, by a rank-one perturbation T = E\d){d\ - that is to say T^(x) = 
(l){0)5{x). The Dirac functional 6 is however bounded on the domain of fco- 
The singular nature of the potential implies that the particle wavehmction will 
be discontinuous at the origin. In particular, such functions will not be in the 
domain of the operator ko- 

We now review the theory of sclf-adjoint extensions of the generator of linear 
translations id on the punctured line M/ {0}. 

2.1 Distributions on discontinuous test functions 

Let AC (I) denote the set of absolutely continuous functions on an open subset 
/ of the real line. We consider some singular functionals on this space. 

Definition 1 The one-sided delta functionals S± on AC (R/ {0}) are given by 

(<5±|(/)) = 0(O±), 

and we introduce the associated functionals: 

(jump at the origin) {j\ = {5+ — (5_| , 
(symmetric delta functional) ((5*| = ^ {S+ -|- 15_| . 

That is, - 0(0") and {6^\(j)) = i</>(0+) + ^0(0"). 

Note that we have 

\S±) = \S.)±l\j). 

Lemma 2 Introducing the form J = |(5+) — (5_|on domain AC (R/ {0}), 
we have the following identity 

J = \j){S4 + \S.){j\. (1) 

Proof. The right hand side is ^\S+ - S-) {5+ + 5- \ + l\S+ + S-) {S+ - 5-\ 
and expanding out gives the result. ■ 

2.2 Distributional First Order Derivatives 

The differential operator ko = id is symmetric onW^'"^ (M) and its closure is 
the generator of translations of wavefunctions along the real axis. This is not 
true when we try restrict to W^'"^ (R/ {0}) due to the jump discontinuity at the 
origin In fact, id is not now symmetric as we have the integration by parts 
formula 

i^pm + {dm = rK- ^ {^\j<^) . 

for (f),il) & W^'"^ (R/ {0}). Instead, let us introduce the distributional derivative 
operator 

iD^id + i\S,){j\. (2) 

A combination of the lemma and the integration-by-parts formula shows that 
iD is symmetric operator on W^'"^ (R/ {0}). Note that is then typically a 
functional on W^'^ (R/ {0}) and not a function. 
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2.3 Albeverio-Kurasov Construction 

Albeverio and Kurasov study rank-one perturbations of self adjoint operators fco 
that are not semibounded, specifically they study formal of the form fcg + Iv') {f\ 
where is a bounded functional on dom{kQ) though the perturbation need not 
be form bounded. The consider the restriction ko of ko to the dense domain 
dom{ko) = {ip £ dom{ko) : {(p\ ip) — 0}, and study the scales of Hilbert spaces 
associated with the adjoint fcg. 

In our case, ko = id with domain W^''^ (M) and, since we consider i5-function 
perturbations, the restricted domain will be taken to be Vb = {'0 G W^''^ (R) : 
^ (0) = 0}. In this kl will be the operator id with domain W^''^ (M/ {0}). 

2.4 Deficiency Subspaces 

Let us introduce the following pair of vectors (j>j. S W^''^ (R/ {0}): 

(/.±(i)-T*e±'l(o,oo)(±t). (3) 

An elementary calculation shows that d \(f>±) — T \4'±) and that {j\4'±) = —i, 
so that iD \<j)±) — T« \4'±) + or 

= {iD±i)-^ \S,). 

More generally iD„ |0_|_) = |0±) + 10- 

It is convenient to realize iD as being the adjoint k^ to the operator kg 
as the restriction of id to the dense domain Vb = {i/' G VF^'^ (R) : ip {0) — O}. 
Note that Vq is a Hilbert space with the Sobolev inner product. The deficiency 
subspaces V± = ker(fcg ± i) are then both one-dimensional and spanned by the 
defect vectors (j)^. respectively: V± — {C(j)j.]. The domain of the adjoint can 
then be written as 

domikl) = Vo ©1,2 V+ ®i,2 V- (4) 

where the three subspaces are mutually orthogonal with respect to the Sobolev 
space inner product. (See, for instance, theorem X.2 [TB].) 
The elements of dom{kQ) can then be represented as 

IV') = iV'o> + c+ + c_|(/)_) 

where ipo ^ and c± = ±i4> (0*) G C, the action of the adjoint is 
kl 1^) = id - ic+ + ic^ . 

We therefore represent iD by fcj and to check consistency, we note that for 
ip G c?om(fcJ) 

tD\i:) = z(}|Vo)+iV'(0+)(-z + |<5.))-z^(0-)(z + 
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For G dom{kl), the associated boundary form [(j)\tp] :— + (^o'/'IV') 

will be [4>\tp] = —i All self-adjoint extensions of ko can then be parame- 

terized by a unimodular complex parameter s. We define fco,s as the restriction 
of kl to the domain 

dom (^ko^s) (v- e dora{kl) : iP (0+) = (O^)} , 

where the functions have an abrupt phase change s as we pass through the 
origin. 



2.5 Singular Perturbations of id 

We return to the formal Hamiltonian k ^ id + E5 which we now interpret as 

k = iD + E\S^) {S^\ 

The Hamiltonian can then be split up into continuous and singular components 
k = kac + A:sing where k^c = id and 

fcsing = j 1^*) 01 + E (5*1 . 

Generally speaking, kli/j) will be a functional for given i/j G dom{kQ) except 
however when fcsing \ip) = and the space of vector for which this vanishes 
defines the domain of k. This requires that {i{j\ + E IV-") = 0, or 

t (0+) - V (0-)] + ^E [V (0+) + ij (0-)] = 0. 
This can be written as the boundary condition 

V'(o-) = iq4|v'(o+). (5) 

The domain of k is therefore the set of ?/> e W^'^ (M/ {0}) satisfying ([5]). Note 

1-^E ^ 
that the phase change s — is the Cayley transform of ^E. 



3 Second Quantization 

We now consider the Hilbert space 

Sj^i)(E,T{Ll (R, dt)) 

where are fixed separable Hilbert spaces and T (•) is the bosonic Fock functor. 
A typical vector in is $ = ($n) where, for n = 0, 1, 2, • • ■ , we have that <i>„ is 
a f) C>5 .S®"-valued function on M" symmetric in all its arguments and such that 

y2~l ,tn)fdti---dtn <00. 

„>o -^R" 

We shall treat the case ^ = C initially for transparency. 
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3.1 Second quantization of iD 

Fock space has the continuous tensor product decomposition F t)2) — 
F (f)i) F ([)2) and we exploit the direct sum decomposition ^ to write 

F (doTO(fcJ)) ^ Fo ®i.2 F+ ®i,2 F_, (6) 

where F^ = F for # = 0, +, — , and the tensor product is with respect to 
the Sobolev inner product. Recall that every -0 G dom(A;J) will have the decom- 
position ^ = ipQ + iip (0^) — iip (0~) (f)_ and can define the usual exponential 
vectors e {tjj) — J2n>o TTT ®i 2 ^ ^-nd this domain define the annihilation fields 
A (0), with G doin(kl), by 

where (<;6|V')i,2 = ('/'olV'o)i,2 + V-* (0+) </> (0+) + V* (0") <^ (0"). 

Now the spaces V± are both one-dimensional and so the Fock spaces F (V±) 
each correspond to the Hilbert space of an independent single mode harmonic 
oscillator and we take the respective annihilator operators to be 

a±^A (±i0±) . (7) 

With this convention, wc have for (p e doni{K'^), A {(j>) = A{(I)q) + cj) (0+) a+ + 
(/) (0~) a_. In particular, since (±i(/)^|'0)j^ 2 ~ (0^)? we have 

a±e (iP) = iP {0^) eW. 
It is convenient to introduce 

a* = - (a+ + a_) 

The second quantization of fcj = iD is then given by Kq — ifo.ac + -K^o.sing 
where 

Ko,^ = dr{W), Ko,sing = idT{\5.){j\) = ml {a+~ a-). (8) 

3.2 Singular Perturbations of K^c = dV (id) 

We now consider the perturbed Hamiltonian K = Kq + T where 

T = E'liala* -t- Eioal + £'01 a* + Eqq (9) 

where the Ec,f3 are operators on t) with (-Eq/j)^ = Epa. The Hamiltonian may 
then be decomposed into continuous and singular parts as 

Ki,c + Eqq + i?oia* + al {ia+ - ia- + Ena^, + Eiq) . 

Let us introduce the following operators on the system space 

l~_jEii 1 H ^ EQo + -Eoilml 1 \ Eiq. 

1 + l + ^En 2 ]^i + iEnj 

(10) 
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Theorem 3 The domain of the self- adjoint extension corresponding to the op- 
erator K = Kq + T is the set of vectors $ S ()®dom(_ftro) satisfying the boundary 
conditions 

a-^ = Sa+<P + L, (11) 

and on this domain we have 

- iHtotai^ = (^-iKac - \l^L -iH^'P- L^Sa+^. (12) 

Proof. Our strategy is basically the same as in the one-particle case - 
we choose the domain of the operator to consist of those vectors $ such that 
(ia^ — ia^ + Ena^, + -Eio) $ = 0. This is equivalent to the boundary condition 
(jlip . For vectors on this domain, we may substitute a_$ for a+$ using the 
boundary condition and this gives the desired result. ■ 

3.3 Multiple Field Modes 

We now consider the more general case where .ft = C". Let {cj : j = 1, - ■ ■ ,n} 
be an orthonormal basis for The one-particle space is now (R, dt) = 

C" (R^ dt) and we replace the decomposition d!]) by T (c"- (g) dom{kl)j = 

r[|"^ 01,2 r^"^ 01,2 rL"^ where = r (C" ® V#). The defect vectors for the 
amphation of ko to C" (R, dt) may now be fixed as Cj (E> 4'± so the 
deficiency indices are now (n,n). 

Proceeding as before, we introduce the independent annihilators 

aj^± = A {±iej (X) 4>±) 

and the second quantization of the ampliation of fcj is Kq — -R'o.ac + -f^o.sing with 

Ko^e^c = dr (id) , Xo,sing = «a],* (aj,+ - . (13) 

(A summation over the range 1, ■ • • , n is implied for repeated Latin indices!) 
Let us also set ao,± = 1. As perturbation wc consider the singular term 

T = aiap. (14) 

(Here we understand repeated Greek indices as implying a sum over 0, 1, • • • , n.) 
Again the {n -\- 1)^ operators Ea(3 are operators on 05 (i)) and we require that 
(Eaji)^ = Epa- It is convenient to assemble them into a matrix E. More 
generally we consider the class of matrices 

X= f ?° ) e»((C®J?)®f)) 

where Xqq G 05 (()) , Xqi, e <B (i), J? ® fi) , e 03 (jl ® f), f)) and Xu G 05 (i^ f)) . 
That is Xqi is the row vector (A"oi, • ■ • , A'on) and X^q is the column vector 
(Xoi, • • • , Xon)^ while Xu = {Xij). 
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Let us also introduce the special matrix 11 projecting onto the subspace ^^t) 
of (C ® J^) (g) f) EE () © (il (g) f)), that is, 



n 




1 



Its coefficients are the Evans-Hudson delta 



I, P G {!,■■■ ,n}, 

0, otherwise. 



Given a self-adjoint *B (f))-valued matrix E as above, we recall some related 
matrices, originally introduced in [13j . starting with the matrix G defined by 
the identity G = -tE-iEHG, 



G + n 

(i + nG) = 
f 



1 
L S 



(Ito) G = + ^En 1 E 

(Model) V 
(Galilean) M 
(Dressing) F 
Note the identity 

G = iEF 

and that we encounter the operators 



L S-1 



^ 1+inG = i+-nE 



f 



(15) 



S 



1 - -Ee 
2 



l + -Eu] G<B(j?«)f)), 



L = -i{l + -Eu] Eeae'B{A®l),l)). 



H 



En 



llmEoe ( 1 + -^;, 



Em. 



(16) 



With these conventions we may write 



— ialll (a+ — a_) , T = a]^Ea*, 



where 



a± 



/ 1 \ 

ai.± 



' a^ = -(a++a_). 
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The Hamiltonian is then 

K ^dr (id) + al {in (a+ - a_) + EaJ 

and the domain is the set of vectors $ such that 11 {i (a+ — a_) + Ea^} <i> = 0. 
This boundary condition ensures that all the ^ terms vanish {j = 1, ■ ■ ■ ,n) 
and may be reformulated as 

n (E - 2m) a+4> + n (E + 2in) a_$ = 

or 

with S = Sjk CS) \ej) {ek\ and L = Lj ® \ej). We have trivially that ao,_ = ao,+ 
and the we may include this in the boundary condition to write 

a+$=Ma_$. (17) 

Using the boundary condition, we may write 

= i^ac$ + ^-Boa {do,p + Map} ap^+^ 

however ^E (1 + M) = E (l + \TIG) = EF = iG. With this we may write the 
action of the Hamiltonian on vectors satisfying the boundary conditions as 

= K^c^ + iGoaap,+^. 

We may summarize our findings using the model matrix V in the following 
theorem. 

Theorem 4 The Hamiltonian associated to K = Kq + T with perturbation 
is defined on the domain T ^dom(kQ)j satisfying the boundary condition 
a+$= Va_<J>. Its action on this domain, and the boundary condition, are given 
by 

~iK^ = (Voo + Vokak,+ - iKo) * = {Vopop (0+) - iK„) $, 

4 Gauge freedom 

A complex number k+ with strictly positive real part will be referred to as a 
complex damping constant. For convenience we shall normalize complex damp- 
ing constants as 

K± — — ± ia, 
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where a is real. For a fixed complex damping k we then define a functional 
C = Ck by 

10 = \5+) + ^ 7 \5,) + |j) . (18) 

We then we have the following local identity generalizing ([T]) 

:^-b)(CI + IC)OI, (19) 

as the \j) {j\ terms cancel. This allows us to construct a more general class of 
self-adjoint extensions of the restriction of id. Let a be the imaginary part of 
the complex damping k+ and define 

iD,^id + i\5,){:)\^a\j) {:j\ = id + i\c) (j| 

with Q as above. It follows that iD„ is likewise a symmetric operator on 
VFi'2(M/{0}). 

The formal Hamiltonian k = id + E5 may be alternatively interpreted as 

K=iD„ + E\C) (CI, 

or fco- = id + i \() {j\ + E |C) (CI- We may follow the same argument as before 
and arrange for the singular component to vanish by imposing the boundary 
conditions: this time the condition ([5]) is modified to -0 (0~) = So-V' (O"*") where 

we now have the phase Sa — ■ 

^ 1 + iK+E 

This can be lifted immediately to the second quantization. For .ft = C, we 
set a = K_a+ + K+a_. We make the corresponding replacements: i^o,sing (c) = 
ia^ (a+ — a_) and T (cr) = Eno^a + Eiqo) + E^ia + Eqo- The boundary con- 
ditions arise from requiring the terms to vanish and after similar algebra 
to before we arrive at the restatement of the first theorem with the modified 
operators 

„ . s 1 — 1 , , i 

= 1^- F^ L{a) = -———E^o, 
1 + m+Eii 1 + iK^Eii 

H{a) = Eoa+lm{Eoi—^^^—Eio\. (20) 

1+lK+Ell J 

Identical expression where obtained for the singular limit of finite time correlated 
Bose field in [12 , theorem 8.1 equation (8.11). 

The multiple field mode version of this is to introduce in place of a^, = 
i (a+ + a_ ) the vector 

a=(il + iZ)a+ + (il-iZ)a_ 
where Z=[[' ^ )isa hermitean matrix with complex entries. We then ob- 



. Zt_ ^ 

tain a straightforward modification of the second theorem with new coefficients 
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G (E, Z) = -i (l + (in + iZ)) E. This general type was first introduced 
in [13]. 

In general, the parameters (Zjk) may be termed "gauge parameters" and by 
choosing a different set of parameters we obtain a different self-adjoint exten- 
sion. This type freedom/ ambiguity is well-known and is ultimately a question 
of fixing the desired physical model and cannot arise from purely mathematical 
arguments alone [16j . 

Acknowledgement 5 The author wishes to thank Matthew James and Oleg 
G. Smolyanov Jar fruitful discussions. He is particularly grateful to Alexei 
lantchenko for introducing him to reference J^. 
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